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We establish an exact partial differential equation to model Kerr comb generation in whispering- 
gallery mode resonators. This equation is a variant of the Lugiato-Lefever equation that includes 
higher-order dispersion and nonlinearity. This spatio-temporal model, whose main variable is the 
total intracavity field, is significantly more suitable than the modal expansion approach for the 
theoretical understanding and the numerical simulation of wide-span combs. It allows us to explore 
pulse formation in which a large number of modes interact cooperatively. This versatile approach 
can be straightforwardly extended to include higher-order dispersion, as well as other phenomena 
like Raman, Brillouin and Rayleigh scattering. We demonstrate for the first time that when the 
dispersion is anomalous, Kerr comb generation can arise as the spectral signature of dissipative 
cavity solitons, leading to wide-span combs with low pumping. 

PACS numbers: 42.62.Eh, 42.65.Hw, 42.65.Sf, 42.65.Tg 



The development of frequency combs - equidistant fre- 
quency lines from a short-pulse laser - revolutionized 
the measurement of frequencies [1] and has opened up a 
host of potential applications in fundamental and applied 
physics, including the measurement of physical constants, 
the detection of earth-like planets, chemical sensing, the 
generation, measurement, and distribution of highly ac- 
curate time, and the generation of low-phase-noise mi- 
crowave radiation pi. Tiisapphire lasers were used as 
the original source of frequency combs, but in the past 
seven years, alternative fiber laser sources have devel- 
oped [2]. Recently, Del'Haye, et al. [5] demonstrated 
that it is possible to use the whispering gallery modes 
in microresonators in combination with the Kerr effect 
to generate an equidistant frequency comb that is also 
referred to as a Kerr comb. Since many applications 
in both fundamental and applied science would benefit 
from the small size, simplicity, robustness, and low power 
consumption of these whispering gallery mode sources, a 
considerable worldwide effort has gone into understand- 
ing and controlling them 0]. In particular, there have 
been several efforts to develop mathematical models of 
these sources [SHE], but all the efforts to date have seri- 
ous drawbacks. 

A complete modal expansion has been derived to de- 
scribe the growth of the Kerr combs from noise [5J [B]. 
This model predicts a cascaded growth in which a pri- 
mary comb is first generated, which then generates a sec- 
ondary comb, and later higher order combs. This model 
is in complete agreement with experiments [5]. However, 
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it is difficult and computationally expensive to use this 
mode expansion beyond the primary comb generation be- 
cause the number of modes that must be kept in a cal- 
culation grows like the third power of time. Moreover, it 
is difficult to study pulse formation in this model, since 
a large number of modes interact cooperatively. Pulse 
formation plays a critical role in comb generation. It is 
desirable to find a spatio-temporal model that would be 
analogous to the Haus modelocking equation (a variant of 
the nonlinear Schrodinger equation) that has been used 
with great success to study modelocked lasers [5]. How- 
ever, this equation is not itself appropriate for whisper- 
ing gallery mode resonators. This equation assumes that 
pulses are isolated, and it does not take into account the 
periodic boundary conditions and the fundamental role 
that continuous wave pumping and dissipation play in 
the evolution [8^. 

In this letter, we show that the Lugiato-Lefever equa- 
tion (LLE) [ro is the appropriate spatio-temporal model 
of whispering gallery mode resonators. We demonstrate 
that the previously-derived modal expansion is equiv- 
alent to a variant of the LLE that includes higher- 
order dispersion and nonlinearity, and the standard LLE 
closely approximates the modal expansion. We then use 
this model to study pulse formation and determine condi- 
tions under which multiple pulses (rolls) and single pulses 
(dissipative solitons) form in the cavity. We show that 
a broadband comb is obtained when a single dissipative 
soliton forms. 

The LLE can be considered a variant of the nonlin- 
ear Schrodinger equation (NLSE) that includes damping, 
driving, and detuning. It has been the subject of inten- 
sive mathematical study pTUl3] , and it has been experi- 
mentally demonstrated that this model can be applied to 
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fiber ring resonators [T3]. Unexpected phenomena such 
has convective instabihty |15j and excitabihty jl6j are 
also known occur in this system. The identification of the 
LLE as the fundamental equation governing the field evo- 
lution in "whispering gallery mode resonators allows us to 
take advantage of this prior work, while at the same time 
suggesting efhcient computational methods for studying 
the field evolution, including Kerr comb generation. 

Our starting point is the modal equations that were de- 
rived in ref. |^. The modal structure of the whispering 
gallery mode resonators is well-understood. It has been 
shown that these resonators can sustain several families 
of eigenmodes that trap light inside the cavity by total 
internal reflection. We are only interested in the fun- 
damental family of modes that can be characterized by 
their toroidal structure. Provided that the polarization is 
fixed, the members of this family are unambiguously de- 
fined by an integer wavcmimber i that characterizes each 
member's angular momentum and can be interpreted as 
the total number of reflections that a photon makes dur- 
ing one round trip in the cavity. 

We denote the eigennumber of the pumped mode as 
Iq. If we restrict ourselves to the spectral neighborhood 
of the eigenfrequencies can be expanded in a Taylor 
series, whose first three elements are 



where 5{x) is the Kronecker delta-function that equals 1 
when a; = and equals zero otherwise. The mode fields 
have been normalized so that \Af\'^ corresponds to the 
photon number in the mode £. The mode bandwidth 
Acj£ = ui/Qq is inversely proportional to the loaded 
quality factor Qo and to the phonon lifetime Tph,^ — 
l/Awf. The four- wave mixing gain is = n2c/iajQ/ngVoj 
where h is Planck's constant, n2 is the Kerr coefficient, 
and Vb is the effective mode volume. The coefficient 
p^m -n p j]2(jicates the mode overlap. The parameter J-q 
denotes the amplitude of the external excitation, while 
flo is the angular frequency of the pump laser, and it is 
assumed to be close to ujq. Typically, resonant pumping 
only occurs when jfio — wo| < Awq- 

The spatiotemporal slowly varying envelope of the to- 
tal field A{0,t) may now be written 

Aie, t) = J2 Mt) exp[i(w, - ujo)t - i{l - , (3) 

where d € [— i", tt] is the azimuthal angle along the cir- 
cumference. From Eq. ([3|, it follows that 
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(^l'^ Equation ([3| independently yields 



where uJe^ = ujq is the eigenfrequency at ^ — £o, 
Ci = dtu/dilg^^^ = AcjpsR is the free-spectral range of 
the resonator or the intermodal angular frequency, and 
C2 = d^uj/di'^\^_^ is the second-order dispersion coeffi- 
cient. The quantity C2 corresponds to ( in refs. |5l [6]. 
In the case of a disk resonator with main radius a, we 
find Ci = c/noa, where c is the velocity of light and no 
is the index of refraction at ujo- This intermodal angular 
frequency is linked to the round-trip period of a photon 
through the resonator as T = 27r/Ci. The second-order 
dispersion ^2 denotes the lowest-order deviation from fre- 
quency equistance of the modes. When (2 = 0, the eigen- 
frequencies are equidistant to lowest order and are sepa- 
rated by Ci. The dispersion is normal when C2 < and 
anomalous when C2 > 0. In fact, C2 is the sum of two 
contributions - the geometrical dispersion (generally nor- 
mal) and the material dispersion (which can be normal or 
anomalous). Explicit expressions are given in refs. [HIS] 
for a spherical resonator. 

We now consider the elements of the modal expansion 
in a range of ^-values in which the expansion of Eq. ([T]) 
is valid. The slowly varying envelopes Ai obey the equa- 
tions E] 
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so that if we restrict ourselves in the degenerate case 
where A^"^"^'' = 1 and Aw^ = Awq, the evolution equa- 
tion Q can be rewritten as 

BA^ 1 

= --Ac.o^-z5o|-4|M+-Ac.o-Foe''^* 
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where cr = fio — is the detuning between the laser 
and cavity resonance frequencies. It is useful to translate 
the frequency of the carrier envelope to remove the ex- 
plicit time-dependence of the driving term by making the 
transformation A^ A exp(icri). It is also useful to trans- 
form the 0-coordinate to remove the group velocity mo- 
tion by making the transformation 0-^0 — Cit mod [27r]. 
When that is done, we find that Eq. ([6]) becomes 
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Finally, this equation can be rewritten in the form of the 
normalized Lugiato-Lefever equation 
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FIG. 1: Multiple-FSR comb with a = 0, /3 = -0.0125 and 
F — 1.71. The intracavity pump is IV'ol^ ~ 1-2, so that it 
is above threshold IV'lth = 1- The parameters correspond to 
Fig. 10b in ref. [6]. (a) Time-domain dynamics that consists of 
19 smooth pulses ("roll" Turing pattern solution). (6) FFT of 
the roll solution, corresponding to an order 19 muItiple-FSR 
comb. 



where the field envelope has been rescaled so that 'ip — 
(2go/Awo)^^^-4* and the time has been rescaled so that 
r = AujQt/2. The dimensionless parameters of this 
normalized equation are the frequency detuning a — 
—2a/AoJo, the dispersion /3 = — 2(^2/Awo, and the ex- 
ternal pump F = (2go/Aajo)^/^^o ■ This LLE with peri- 
odic boundary conditions is the exact counterpart of the 
modal expansion as long as higher-order dispersion and 
the variation of Aoje and A^'"^"^'' can be neglected. Since 
the higher-order corrections can be calculated, it is al- 
ways possible to check the validity of the LLE, and it will 
remain valid until the comb is nearly octave-spanning. As 
a consequence, all the conclusions that were previously 
obtained using the modal expansion and confirmed ex- 
perimentally can also be obtained by solving the LLE. 
These two twin-models are however useful in different 
and complementary ways. On the one hand, the modal 
expansion must be used to determine threshold phenom- 
ena when a small number of modes are involved. On the 
other hand, the LLE is appropriate to use when many 
hundreds or thousands of modes interact since it does 
not refer to the individual modes. In particular, as we 
will show later, it is useful in the study of Kerr combs or 
pulse (soliton) growth and propagation in which a large 
number of modes interact cooperatively. 

We now discuss the meaning and order of magnitude 
of the variables and parameters that appear in the di- 
mensionless LLE. The detailed analysis of ref. [B] showed 
that the threshold value for the Kerr comb is given by 
\A\^Y^ = Acjo/2.go, which exactly corresponds to j-^j^j^ = 1. 
Hence, we would expect ijj to be on the order of 1 in our 
theoretical analysis and numerical simulations. Note that 
the stationary intra-cavity field -00 before comb genera- 
tion can be obtained by setting all the derivatives to zero 
in Eq. (|8]). The dimensionless time can be rewritten as 
i/2rph, where we recall that Tph is the photon lifetime, 
which is typically a few /is in ultra-high-Q whispering 
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FIG. 2: Temporal cavity soliton with a ^ 1.7, P ^ -0.002 
and F = 1.22. The intracavity pump is |i/'o|^ = 0.95, so that 
it is below threshold, (a) Cavity soHton formation, (b) FFT of 
the soliton, corresponding to a wide-span, low-threshold Kerr 
frequency comb. 



gallery mode resonators. Using the external field's phase 
as our reference, the amplitude F is real and positive and 
F^ is proportional to the external pump power. The real 
parameter a equals the ratio of the detuning to half the 
linewidth, so that we expect to be off resonance when 
\a\ > 1. This parameter is easily tuned in experiments. 
The parameter f3 equals the ratio of the walkoff due to 
second-order dispersion to the half-linewidth. This pa- 
rameter is negative for anomalous dispersion and positive 
for normal dispersion. It is a relatively small parameter. 
In refs. O |6] for example, it was found that /3 ^ —0.01 
was sufficient to generate combs. In practice, the disper- 
sion magnitude \f3\ should not be too large (~ 1), as it 
would mean that the deviation from equidistance of the 
mode frequencies would be too strong for the nonlinearity 
to compensate, and the emergence of a wide-span Kerr 
comb would be suppressed. 

Numerical simulations were performed using the split- 
step Fourier algorithm, which is commonly used in sim- 
ulations of the 1-D generalized NLSE. We note that this 
simulation method inherently assumes periodic boundary 
conditions because it is based on the fast Fourier trans- 
form (FFT). For studies of pulses in optical fibers - in 
which the pulse duration can often be many orders of 
magnitude smaller than the separation between pulses - 
this feature is often an annoyance that leads to boundary- 
induced artifacts. In our case, the physical boundary con- 
ditions are periodic. The split-step Fourier algorithm is 
therefore a remarkably cost-effective computational tool, 
enabling the simulation of Kerr comb dynamics with a 
laptop computer in a few minutes regardless of its spec- 
tral span, as opposed to a few days with the modal ex- 
pansion for wide-span combs. 

For the numerical simulation of the LLE, Eq. ([s]), we 
have considered a calcium fluoride resonator with a ra- 
dius a — 2.5 mm. The polar eigennumber of the TE- 
polarized pump mode is £o = 14350, corresponding to a 
wavelength of 1560.5 nm and ujq ^ 2?: x 192.24 THz in 
vacuum. The index of refraction is uq — 1.43. The res- 
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onator is critically coupled with a loaded quality factor 
Qo = 3 X 10^, corresponding to a central mode band- 
width of Auq = wq/Q — 27r X 64 kHz. The free spectral 
range is equal to AwpsR = Ci = 27r x 13.36 GHz. 

In Fig.[T] we present simulation results that correspond 
to those that have already been obtained from the modal 
expansion, namely the multiple-FSR solution. It can be 
seen that the solution corresponds to the formation of 
"rolls" in the time domain. We take advantage of the 
LLE formulation to show particular solutions that would 
be difficult to observe in a strictly modal study. Fig- 
ure [2^ shows the formation of a single-peaked cavity soli- 
ton. This dissipative localized structure is a sub-critical 
Turing pattern. The figure corresponds to a single pulse 
of width AT ~ 500 fs, circulating inside the cavity with 
a round-trip time T = 27r/Aa;FSR — 75 ps. In the fre- 
quency domain, this pulse consists of several tens of mod- 
elocked whispering gallery modes, whose frequencies have 
been nonlinear ly shifted so that they are equidistant, as 
shown in Fig. [2]d. Because this cavity soliton is subcriti- 
cal, it emerges abruptly and at a low pump power. The 
solitonic Kerr combs do not grow as the pump power in- 
creases; instead, they are destroyed. By contrast, super- 
critical combs like the one shown in Fig. [l] have spectral 
components that grow in number and power as the exter- 
nal pump power grows. However, the dissipative cavity 
solitons are robust, and their full width at half maximum 
decreases with the dispersion parameter \/3\. 

For the simulation of wide-span combs, the LLE can 
be generalized to take into account higher-order dis- 
persion {(n = d^^/d£"\e^i„ 7^ 0, with n > 3), non- 
degeneracy of the bandwidths (Aw^ ^ Awq), and spatial 
non-degeneracy of the eigenmodes (A^"^"^" ^ 1). In- 
deed, these higher-order effects have to be taken into con- 
sideration for the simulation of octave-spanning combs. 
If we account for third-order dispersion and the first- 
order variation of A^"^"^'' , Eq m\ becomes 
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where r]( = dA/d£ for £rn = = (-p = i = io, with 
r]e^, r]i^, and rje^ being defined analogously. The details 
of the derivation are analogous to similar derivations of 
the NLSE with higher-order corrections |T7]. When the 
bandwidth of the resonator field approaches the carrier 
frequency, the Taylor expansion in Eq. ([I]) is no longer 
valid, and one must use the full functional dependence 
of uJi — ujo on £, by analogy to what is done in studies of 
supercontinuum generation. 



In conclusion, we have demonstrated that the Kerr 
comb evolution in whispering gallery mode resonators 
can be modeled using the Lugiato-Lefever equation and 
its extensions. We have shown that low-threshold, wide- 
span combs can emerge as dissipative cavity solitons. 
With different parameters, rolls appear, corre- sponding 
to narrower-band combs. So, when the goal is to opti- 
mize the bandwidth a correct parameter choice is criti- 
cal. Additionally, we demonstrated that the LLE can be 
extended to incorporate higher-order dispersion and non- 
linearity, and we expect that it can be extended without 
too much difficulty to include Rayleigh, Brillouin, and 
Raman scattering [TH]. A key advantage of the spatio- 
temporal model that we have developed is that the LLE 
has already been the subject of extensive mathematical 
study, and it should be possible to take advantage of this 
earlier work to shed additional light on the dynamical 
properties of Kerr combs. We expect that a better under- 
standing of Kerr comb generation in whispering gallery 
resonators will be the result, leading to new resonator 
designs that can produce octave-spanning combs. 



Y. K. C acknowledges financial support from the Eu- 
ropean Research Council through the project NextPhase 
(ERC StG 278616). C. R. M acknowledges hospitality 
and support from the FEMTO-ST Institute within the 
framework of the LabEx Action. C. R. M. also acknowl- 
edges useful discussions with T. Sylvestre. 



[1] J. Ye, and J. Hall, Scientific American 294, 74-81 (2008). 
[2] S. A. Diddams, J. Opt. Soc. Am. 5 27, B51-B62 (2010). 
[3] P. Del'Haye, A. Schliesser, A. Arcizet, R. Holzwarth, and 

T. J. Kippenberg, Nature 450, 1214-1217 (2007). 
[4] T. J. Kippenberg, R. Holzwarth, and S. A. Diddams, 

Science 322, 555-559 (2011). 
[5] Y. K. Chembo, D. V. Strekalov, and N. Yu, Phys. Rev. 

Lett. 104, 103902 (2010). 
[6] Y. K. Chembo and N. Yu, Phys. Rev. A 82, 033801 

(2010). 

[7] I. H. Agha, Y. Okawachi, and A. L. Gaeta, Opt. Express 



17, 16209-16215 (2009). 
[8] A. B. Matsko, A. A. Savchenko, W. Liang, V. S. Ilchenko, 
D. Seidel, and L. Maleki, Opt. Lett. 36, 2845-2847 
(2011). 

[9] H. A. Haus, IEEE J. Sel. Top. Quantum Electron. 6, 
1173-1185 (2000). 
[10] L. A. Lugiato and R. Lefever, Phys. Rev. Lett. 58, 2209- 
2211 (1987). 

[11] A. J. Scroggie, W. J. Firth, G. S. McDonald, M. Tlidi, R. 
Lefever, and L. A. Lugiato, Chaos, Solitons, and Fractals 
4, 1323-1354 (1994). 



5 



[12] T. Miyaji, I. Olmishi, and Y. Tsutsumi, Physica D 239, 

2066-2083 (2010). 
[13] G. Kozyreff, Physica D 241, 936-946 (2012). 
[14] F. Leo, S. Coen, P. Kockaert, S.-P. Gozra, P. Emplit, and 

M. Haelterman, Nat. Phot. 4, 471-476 (2010). 
[15] M. Santagiustina, P. Colet, M. San Miguel, and D. Wal- 



gracf, Phys. Rev. Lett. 79, 3633-3636 (1997). 
[16] D. Gomila, M. A. Matias, and P. Colet, Phys. Rev. Lett. 

94, 063905 (2005). 
[17] G. P. Agrawal, Nonlinear Fiber Optics (Academic, 2001). 
[18] R. W. Boyd Nonlinear Optics (Academic, 2008). 



